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Abstract 

Space charge effects can be very important for the dynamics of intense particle beams, as they 
repeatedly pass through nonlinear focusing elements, aiming to maximize the beam's luminosity 
properties in the storage rings of a high energy accelerator. In the case of hadron beams, whose 
charge distribution can be considered as "frozen" within a cylindrical core of small radius com- 
pared to the beam's dynamical aperture, analytical formulas have been recently derived l| for the 
contribution of space charges within first order Hamiltonian perturbation theory. These formulas 
involve distribution functions which, in general, do not lead to expressions that can be evaluated in 
closed form. In this paper, we apply this theory to an example of a charge distribution, whose effect 
on the dynamics can be derived explicitly and in closed form, both in the case of 2-dimensional as 
well as 4-dimensional mapping models of hadron beams. We find that, even for very small values 
of the "perveance" (strength of the space charge effect) the long term stability of the dynamics 
changes considerably. In the flat beam case, the outer invariant "tori" surrounding the origin dis- 
appear, decreasing the size of the beam's dynamical aperture, while beyond a certain threshold the 
beam is almost entirely lost. Analogous results in mapping models of beams with 2-dimensional 
cross section demonstrate that in that case also, even for weak tune depressions, orbital diffusion 
is enhanced and many particles whose motion was bounded now escape to infinity, indicating that 
space charges can impose significant limitations on the beam's luminosity. 
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I. INTRODUCTION 



One of the fundamental problems concerning the dynamics of particle beams in the stor- 
age rings of high energy accelerators is the determination of the beam's dynamical aperture, 
i.e. the maximal domain containing the particles closest to their ideal circular path and for 
the longest possible time. For example, "flat" hadron beams (experiencing largely horizontal 
betatron oscillations) can be described by 2-dimensional (2D) area-preserving maps, where 
the existence of invariant curves around the ideal stable path at the origin, guarantees the 
stability of the beam's dynamics for infinitely long times . The reason for this is that 
the chaotic motion between these invariant curves is always bounded and the particles never 
escape to infinity. Other important phenomena that arise in this context is the presence of a 
major resonance in the form of a chain of islands through which the beam may be collected, 
or the existence of an outer invariant curve surrounding these islands, which serves as a 
boundary of the motion and thus allows an estimate of the beam's dynamical aperture. 

On the other hand, hadron beams with a 2 - dimensional cross section require the use 
of 4-dimensional (4D) symplectic mappings for the study of their dynamics J, [J O, B, Q . 
In fact, if longitudinal (or synchrotron) oscillations also need to be included the mappings 
become 6-dimensional. In such cases, the problems of particle loss are severe, as chaotic 
regions around different resonances are connected, providing a network of paths through 
which particles can move away from the origin, and eventually escape from the beam after 
sufficiently long times. 

In this Letter, we add to these issues the presence of space charges within a core radius 
r c , which is small compared to the beam's dynamical aperture. In other words, we will 
assume that our proton (or antiproton) beam is intense enough so that the effect of a charge 
distribution concentrated within this core radius cannot be neglected. Furthermore, we will 
consider this distribution as cylindrically symmetric and "frozen" (i.e. time independent, so 
that it may be self consistent with the linear lattice) and study the dynamics of a hadron 
beam as it passes repeatedly through nonlinear magnetic focusing elements of the FODO 
cell type. This system has been studied extensively in the absence of space charge effects in 
and the question we raise now is whether its dynamics can be seriously 
affected if space charges are also taken into consideration. 

Space charge presents a fundamental limitation to high intensity circular accelerators. 
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Its effects are especially important in the latest designs of high-intensity proton rings, which 
require beam losses much smaller than presently achieved in existing facilities. It is there- 
fore necessary to estimate the major space charge effects which could lead to emittance 
growth and associated beam loss |9(. The interplay between nonlinear effects, typical of 
single-particle dynamics, and space charge, typical of multi-particle dynamics induced by 
the Coulomb interaction, represents a difficult challenge. To understand better these phe- 
nomena, an intense experimental campaign was launched at the CERN Proton Synchrotron 
[loj | . It is very important, therefore, to develop analytical techniques which could be utilized 
in order to study and localize the associated web of resonances (se e e.g . to obtain an 
analytical estimation of the dynamic aperture, as suggested e.g. in 12|. 

In a recent paper, Benedetti and Turchetti [| used first order canonical perturbation 
theory to obtain analytical expressions for the jump in the position and momenta due to the 
multipolar kicks in such maps, showing that the space charges effectively modulate the tune 
at every passage of the particle through a nonlinear element of the lattice. In particular, 
they derived the new position and momentum coordinates after the nth passage through a 
FODO cell in the thin lens approximation, as the iterates of the 2D map 







\ Pn+1 J 





cosf2(J) — sinf2(J) 
sinf2(J) cosfi(J) 




x " , n = 0,l,2,..., (1) 



where 

J=^p, (2) 
for k = 3, i.e. in the case of sextupole nonlinearities and 

The variables P, X and the parameters entering in the above expressions are related to the 
corresponding ones of the original Hamiltonian 



Z=l,2,.. \ c / 



(4) 



by the formulas 



P = p/u 2 , X = x/u , R c = r c u~ 3/2 , u 2 = u 2 - \, (5) 

' c 
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where u is the depressed phase advance at the center of the charge distribution, p = dx/ds, 
s is the coordinate along the ideal circular orbit, g2{t) is given by 



rt / r 2 



g 2 {t) = u 1 g 1 (u)du , gx{t) = / g(u)du , g ( — ) = -nr 2 c p(r) (6) 
Jo Jo 



p(r) satisfies J °° np(r)dr 2 = 1, 7rr 2 p(0) = 1 and Qp(r) represents the radial charge density. 
Note that if m, q, and Vo denote the mass, charge and velocity of our non-relativistic particles, 
the "perveance" parameter £ = 2qQjmv^ determines the tune depression in (|5]L which must 
be small for the above analysis to be valid lj. 

The stage is now set for the investigation of space charge effects on the dynamics. How- 
ever, the space advances Q(J) needed in ([I]) at every iteration depend on integrals of the 
distribution function g(u) that may not be available analytically. To overcome this difficulty, 
we choose in section II a particular form of g(u) for which these integrals can be explicitly 
carried out not only for 2D maps of the "flat" beam case, but also for 4D maps describing 
vertical as well as horizontal deflections of the beam's particles. 

Thus, in section III we perform numerical experiments to examine the influence of space 
charges on the dynamics and find indeed that even for small perveance values the long 
term stability of the beam is significantly affected. In particular, as £ grows (or the tune 
depression uj/ujo decreases from 1), perturbations of 2D as well as 4D maps show that the 
outer invariant "surfaces" surrounding the ideal circular path at the origin disappear and 
the beam's dynamical aperture is seriously limited. Only the major unperturbed stable reso- 
nances survive, with their "boundaries" clearly diminished by the presence of new resonances 
due to space charge effects. In our 2D mapping model, a threshold value of £ (or uj/ujq) was 
found, beyond which the beam is practically destroyed. The paper ends by describing our 
concluding remarks and work in progress in section IV. 

Space charge effects on beam stability became a relevant topic during the years when 
construction of medium-low energy hig h currents accelerators, such as SNS and the design of 



the FAIR rings at GSI were started 
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13, 



14 1 (see also many articles in the SNS Accelerator 



Physics Proceedings of the last few years). The role of collective effects and resonances has 
attracted considerable attention, since they can cause significant beam quality deterioration 
and losses Another relevant issue is the coupling with the longitudinal motion which 
modulates the transverse tune and induces losses by resonance crossing as shown by recent 
experiments jiol ]. 



High intensity rings, where the bunches can circulate over one million turns, require a 
careful analysis of the long term stability of the beam. Since the commonly used codes 
require large CPU times and exhibit an emittance growth due to numerical noise, they 
are not suited for long term dynamic aperture studies and the use of faster methods is 
necessary ^j]. The method proposed in jl| allows us to introduce space charge effects in 
one single evaluation of the map, when a thin sextupole or octupole is present, just as one 
does in the absence of space charge, and is thus especially well suited for dynamical aperture 
calculations. 



II. EXACT RESULTS FOR A SPECIFIC CHARGE DISTRIBUTION 



A. The One - dimensional Beam 



Let us choose for our space charge distribution function in the 1-dimensional case the 
form 

9 ( 1?) = WM+W ' (7) 

The generalization to 2 dimensions is evident by replacing X 2 by X 2 + Y 2 in (JJJ). Observe 
that this function satisfies the requirements that g(0) = 1, gi(t) oc t as t — > and gi{t) — > 1 
as t — > 0, (using ©), as expected from the theory 

Evaluating now by elementary manipulations the integral in (J3J), using (J7J) and J6J, we 
find that it is given by the closed form expression 

2n /2Jcos 2 6\ 2?r 

Thus, the phase advance at every iteration becomes 

where J is given by (J2J). Note that, in the limit J — > 0, eq. Q implies that Q —>■ uj as 
expected. In fact, expanding the square root in that limit we find 

o/tn ooi-u 2 /3 J 5 J 2 \ 

from which we can estimate the deviation of Q from the depressed tune uo near the origin. 
In section III below we pick an ujq such that for £ = we have a major resonance and a 



relatively large dynamical aperture in the X n ,P n plane, select r c small compared with this 
aperture and vary £ to study the space charge effect on the dynamics. 



B. The Two - dimensional Beam 



Let us now observe that in two space dimensions the original Hamiltonian of the system, 
(HJ), becomes 

£ fx 2 + y 2 ' 



_ p 2 , x 2 p 2 2 y 2 
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where sextupole nonlinearities involve, of course, both x and y variables. Since there are 
now two tune depressions 



2 n 1/2 ( 2 ^ 1/2 



(12) 



after transforming to new variables X = xu^ 2 , P\ = p\OJ x 1 ^ 2 and Y = yu^ 2 , P 2 = p 2 co> 2 ^ 2 
defined by 

X = (2J 1 ) 1/2 cos6 1 , P 1 = -(2J 1 ) 1/2 sm6 1 , Y = (2J 2 ) 1/2 cos0 2 , P 2 = -(2J 2 ) 1/2 sin 2 , 

(13) 

differentiating the Hamiltonian with respect to J\ and J 2 and integrating over d\ and 6 2 , we 
find the two tune depressions 

^ ' 2 *2cxx?9 1 



2w! V' (2vr) 2 io 7o 4 + 1 
and O2, with 1 replaced by 2 in ((H)) , while A is defined by 



-<f0id0 : 



_ 2Jicos 2 0i 2J 2 cos 2 2 

A — 1 



r 2 



r 2 



(14) 



(15) 



i, 1/2 1/2 

where r\ = r c uj{ , r 2 = t c uj 2 . 



Observe that we have used for the g\ function under the integral sign (see (j3J)), the 
expression g\(A) = A/(l + A), derived from our simple choice of the distribution function 
(JZJ) using (0). 

The above flj. and Q 2 are to be used in the iterations of the 4D mapping: 

{ X(n + 1)\ / cos O x -sin ^1 \ 

P x (n+1) sinfii cosfii 

K(n + 1) cosQ 2 -smfl 2 

\P 2 (n + l) J \ sinfi 2 cosfi 2 / 



(16) 



/ X(n) \ 

P 1 (n)+X 2 (n)-Y 2 (n) 
Y(n) 

\ P 2 (n)-2X{n)Y(n) ) 

whose dynamics has already been studied extensively in jj, 0, |fj |2t |8j in the absence of space 
charge effects, i.e for u% = uqx and u 2 = ujq 2 . In these papers, it was observed that for the 
tune values q x = 0.61903, and q y = 0.4152, in 



W01 = 2Tcq x , u 02 = 2nq y , 



(17) 



a large dynamical aperture is achieved, with interesting chains of resonant "tori" surrounding 
the origin. In section III we will study what happens to these structures when £ > (i.e 
uj\ < loqi, uj 2 < uj 02 ) and space charge effects are taken into account. Before doing this, 
however, it is necessary to describe how the integrals in (|T%j) are to be evaluated: Let us 
first perform the integration over 9 2 , writing 



2tt i-2ti 

Jo 



o 2/1 /*27r 

1 d9 1 d9 2 = 2 dBx cos 2 ejtfx 



where 



I (9i 



2- 



d6o 



o 2A X cos 2 6 2 + B 1 



A\ = -| , Bx 



2J X cos 2 9x 



(18) 



(19) 



The integral f)19|) can be evaluated as before with elementary functions to yield 







2tt 



[Bx(Bx+2Ax)} 1 / 2 



2tt 



1 + 



2Jx cos 2 9x 



2Jx cos 2 Bx 2J 2 
1 + — ^ + " 



2 ( 1 1 2 1 2 

r l / V r l r 2 



-1/2 



Inserting now expression (}2T)|) into the integral (|18)L changing integration variable to 
we easily arrive, after some simple manipulations, to the expression 

cosd) + 1 



^1 



2vrr 2 ^ 



where 



Jx J Y [(cos^ + C 1 )(cos(f) + D 1 ))] 1 / 2 ' 

cx=i + i, Dl =i + iu + ^ 



Jx 



Jx 



(20) 

(21) 
(22) 



We finally make the substitution u = tan(0/2) and rewrite the above integral in the form 
167rr 2 f°° du 



Kx 



Ji Jo {u 2 + l){ [{Cx - l)u 2 + 1 + Cx] [{Dx - l)u 2 + 1 + D x \ } 



1/2 



(23) 



7 



This is clearly not an elementary integral. Notice, however, that all terms in the denominator 
of (J23|) are positive and as u — > oo the integrand vanishes as u~ 4 . It is, therefore, expected to 
converge very rapidly and may be computed, at every iteration of the map, using standard 
routines. For practical purposes, however, in section III below, we prefer to compute instead 
its equivalent form (|21|). Of course, as explained above, a similar integral, K 2 , also needs to 
be computed (with 1 — > 2 in (JTSJ), (J2U) and ((22$), whence Q± and Sl 2 are found and the next 
iteration of the 4D map (jl6j) can be evaluated. 



III. NUMERICAL RESULTS 



Let us now turn to some practical applications of the above theory to specific problems 
concerning the stability of hadron beams passing through FODO cell magnetic focusing 
elements and experiencing sextupole nonlinearities, as described in . First, we shall 

consider the flat beam case (JTJ), for the specific tune value q x = 0.21 corresponding to 
frequency u>o = 2nq x = 1.32, exhibiting, in the absence of space charge perturbations, the 
phase space picture shown in Figure ^a) below. As we see in this figure, the region of 
bounded particle motion extends to a radius of about 0.54 units from the origin. There are 
also 5 islands of a major resonance surrounded by invariant curves (or ID "tori"), whose 
outermost boundary delimits the so - called dynamical aperture of the beam. Outside that 
domain there are chains of smaller islands (representing higher order resonances) and chaotic 
regions through which particles eventually escape to infinity. This escape occurs, of course, 
at different speeds due to the well - known phenomenon of "stickiness" , depending on how 
close the orbits are to the invariant curves surrounding the islands. 

Let us now consider a space charge distribution of the form ((7j) with a "frozen core" of 
radius r c = 0.1, which is small compared with the radius of the beam's dynamical aperture. 
Our purpose is to vary the value of the preveance £ > 0, cf. (J1J), starting from £ = 0, to 
estimate the effects of space charge on the dynamics. 

Setting £ = 0.001 {u/ojq ~ 0.97), for example, which is quite small compared with 
r\ = 0.01, we observe in Figure^b) that the picture has significantly changed. In particular, 
the 3-dimensional character of the dynamics (due to the variation of the space advance £l(J)) 
has turned the invariant curves into "surfaces" and has led to the dissolution of the outer 
ones surrounding the five major islands. Furthermore, most of the chains of smaller islands 
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FIG. 1: (a) The phase space picture of a flat beam, resulting from the iterations of the map (^Q), 
with tune q x = 0.21 and frequency ujq = 1.32, in the absence of space charges, (b) The same 
picture with £ = 0.001 (or uj/ujo ~ 0.9426) in the Hamiltonian Notice the dissolution of the 
outer invariant curves and the disappearance of many of the smaller islands of (a) away from the 
origin, leading to a significant decrease of the beam's dynamical aperture. Each initial condition 
is followed for N = 10 4 iterations. 

of Figure [T^a) have disappeared due to the new resonances caused by the presence of space 
charges. To see how all this affects the dynamical aperture of the beam as a function of the 
tune depression uj/uoq we now perform the following experiment: 

Forming a grid of initial conditions of step size Ax = Ap = 0.01 within a square [—1, 1] x 
[—1,1] about the origin (X n = P n = 0), we use ((TJ) to iterate for different £ > (or u/u < 1) 
all points falling within circular rings of width Ar = 0.01 for iV = 10 5 and N = 10 6 iterations 
and plot in Figure El the last r = (X^ + P 2 ) 1 ^ 2 = r esc value, at which at least one orbit was 
found to escape from the next outer ring. The results demonstrate that already at £ = 0.001 
(cg>/cg> — 0.97) our estimate of the dynamical aperture r esc has fallen from 0.54 to about 
0.43. In fact, it remains close to that value (rising somewhat to about 0.5) until £ ~ .006 
(lj /ljq ~ 0.81), where it experiences a sudden drop to r esc ~ 0.03, and the beam is effectively 
destroyed. Of course, once one orbit escapes most of them quickly follow within the next 
one or two circular rings. Note also that increasing the number of iterations from iV = 10 5 

9 



to N = 10 6 does not appreciably change the results, until the sudden drop occurs. 



0.6 
0.5 
0.4 

„ 0.3 
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0.000 0.002 0.004 0.006 0.008 

FIG. 2: Plot of dynamical aperture "radius" r esc estimates of the flat beam case, resulting from 
the iterations of the map (fTjl. with frequency ujq = 1.32, in the presence of space charges, i.e for 
increasing values of £ in (j3J). Note the quick decrease of the aperture value by about 20 percent, 
which remains nearly the same until £ ~ 0.006 {oj/ojq ~ 0.81), when a sudden drop occurs, probably 
due to the appearance of a major new resonance. 

This dramatic change at uj/ujq ~ 0.81 is most probably due to the presence of a major 
new resonance caused by the space charge perturbation. It may be an important effect, 
however, since it occurs at a £ value which is still smaller than the r\ = 0.01 radius of the 
charge core. Of course, long before this happens, already at £ ~ 0.0002 (or uj/ujq ~ 0.994), 
the effective aperture of the beam has been significantly reduced by about 20 percent from 
its value at £ = 0. 

Finally, let us turn to the case of the 4D map (fTrjj) . describing the more realistic case 
of a beam whose particles experience horizontal as well as vertical displacements from the 
ideal path, see (JTTJ). For comparison purposes, we choose the same parameter values as in 
our earlier papers 0, Q, 0], i.e horizontal and vertical tunes q x = 0.61903, q y = 0.4152 
respectively, yielding the unperturbed frequencies (fTTj) used in the mapping equations. In 
Figure a), we iterate many initial conditions X(0), Pi(0), F(0), P 2 (0) around the origin 
and plot on a X(n),Pi(n) projection a picture of the dynamics, for |F(n)| < 0.04, in the 
absence of space charges, i.e. with oji = Uoi, % = 1,2. 
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FIG. 3: (a) An X, Pi projection of a beam with x and y deflections, resulting from the iterations 
of the map t|lt>|). with tunes q x = 0.61903,^ = 0.4152, in the absence of space charges. Note the 
presence of a large region of invariant 2D "tori" about the origin and the 6 "islands" of a stable 
resonance, shown here in a cut of the 4D space, with |y(n)| < 0.04. (b) The same picture with 
£ = 0.0002 in Hamiltonian (fTTj) (or uj\/ujq\ ~ 0.9993, LO2/UJ02 — 0.9985). Each initial condition is 
followed for N = 10 4 iterations. Notice the dissolution of outer invariant curves surrounding the 
origin and the disappearance of the chain of islands of (a) , leading to a significant decrease of the 
beam's dynamical aperture. 



Note the region of invariant tori and a chain of 6 "islands" corresponding to a stable reso- 
nance. Strictly speaking, the motion between these tori need not be bounded as 2D surfaces 
do not separate 4D space and Arno'ld diffusion phenomena could in principle carry orbits 
far away from the origin. However, as has been explicitly shown for this model in 0, 0, Q| > 
such phenomena are extremely slow and hence a domain with radius of the order of 0.5 can 
be effectively considered as the dynamical aperture of the beam. Repeating this experiment 
in the presence of space charges, i.e. with £ = 0.0002 (or ui/u i ~ 0.9993, ^2/^02 — 0.9985) 
in (jllj) . we observe in Figure Efb) that the outer invariant curves (together with the islands) 
have been destroyed and the dynamical aperture of the beam has been significantly reduced. 

Studying this reduction as a function of £, we proceed to choose initial conditions from 
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a grid of step size 0.05, extending from -0.65 to 0.65 in all 4 directions about the origin, 
represented by X(0), Pi(0), Y(0), ^2(0). Iterating the resulting orbits from points within 
spherical shells of width Ar = 0.01, we plot in Figure HI for each £, the r esc value of 
the inner radius of the shell from which at least one orbit escapes to infinity. Our results 
demonstrate that the beam's dynamical aperture steadily decreases as £ grows. At £ = 0.006 
(or cui/c^oi — 0.98, UO2/0002 — 0.955) its radius has fallen by more than 50 percent from its 
original value, while at higher perveance values the approximation £ << r\ no longer applies. 
In fact, it is worth noting that the size of the dynamical aperture falls drastically even for 
small values of £, as our calculations with N = 10 5 iterations show. For example even at 
£ = 0.0002 (or uji/ujqi ~ 0.9993, (^2/^02 — 0.9985) our estimate of the dynamical aperture 
has dropped from 0.54 to 0.37. 

0.6 
0.5 
0.4 
„ 0.3 
0.2 
0.1 
0.0 

0.000 0.002 0.004 0.006 0.008 0.010 

k 

FIG. 4: Plot of dynamical aperture estimates r esc of a beam with 2-dimensional cross-section, 
resulting from N = 10 4 and N = 10 5 iterations of the 4D map Q16|). with tunes q x = 0.61903, q y = 
0.4152 for increasing values of £ in pi|). 




IV. CONCLUSIONS 

High intensity effects have long been studied in connection with the so called beam-beam 
interaction and were a relevant topic in the design of many hadron colliders like ISABELLE 
and the SSC (see articles in [l^, [3, 13 )■ However, the effects of high currents on the 
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beam stability have become especially crucial only in recent times, when the design and 
construction of medium energy high current accelerators has started. 

We have reported in this Letter our results on the possible importance of space charge 
effects to the global stability of intense hadron beams, experiencing the sextupole nonlin- 
earities of an array magnetic focusing elements, through which the particles pass iV = 10 5-6 
times in a typical "medium term" experiment of intense beam dynamics. We have used a 
recently developed analytical approach [| to model the space charges by a "frozen core" 
distribution, valid to first order in canonical perturbation theory. By proposing a simple ex- 
ample of such a distribution, which leads to explicit and convenient formulas, we have been 
able to carry out detailed numerical investigations on perturbations of 2D and 4D mapping 
models, describing the dynamics of flat (horizontal) and elliptic beams (with horizontal and 
vertical displacements) respectively. 

These charge distributions are in effect periodic modulations of the tunes (and space 
advance frequencies) of the motion and are therefore expected to introduce new resonances, 
raising the phase space dimensionality of the dynamics. Thus, outer invariant tori of the 
unperturbed case start to disappear and "island" chains of higher order resonances far from 
the origin eventually drift away, leading to a significant decrease of the region of bounded be- 
tatron oscillations of the particles about their ideal path (i.e. the beam's dynamical aperture, 
or luminosity). 

In our experiments, we have been able to measure this reduction of the beam's dynamical 
aperture, for several small values of the perveance parameter £, representing the strength of 
the space charge distribution. We found that, within the range of validity of our approxima- 
tions, the domain of bounded orbits decreases by a significant percentage and hence space 
charge effects should be taken into consideration as they can be important for the long term 
survival of the beam. In the flat beam case, we observed a near total loss of the beam at 
some £ value, which is most likely caused by the onset of a major new resonance introduced 
by the space charge modulations. On the other hand, in the more general case of a beam 
with 2- dimensional cross section modelled by a 4- dimensional map, we also discovered a 
sudden drop in the dynamical aperture, occurring already at very small tune depressions. 

We, therefore, believe that space charges are important enough to merit further investi- 
gation in mapping models of intense proton beams jscj. The occurrence of new low order 
resonances poses, of course, a major threat to the dynamics, if the perveance parameter is 
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big enough. However, even at small values of this parameter, weak (Arnol'd) diffusion effects 
and the slow drift of high order resonances, may significantly alter the long term picture 
of the motion, after a sufficiently great number of iterations. It would also be useful to 
compare the one turn map with the full integration of the space charge effect over one beam 
revolution to appreciate the validity limits of our approximation. Indeed, since the high 
computation efficiency of the one turn map is a key issue of this approach, an estimate of 
the errors in some reference cases would contribute additional useful information in realistic 
applications. 
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